Quantum phase transition in the delocalized regime of the spin-boson model 
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The existence of the delocalized-localized quantum phase transition (QPT) in the ohmic spin- 
boson model has been commonly recognized. While the physics in the localized regime is relatively 
simple, the delocalized regime shows many interesting behaviors. Here we reveal that in this regime 
there exists a novel QPT: namely, from a phase without a bound state to a phase with a bound 
state, which leads to completely different dynamical behaviors in these two phases. If the reservoir is 
initially in the displaced vacuum state (i.e., the coherent state), the spin dynamics exhibits lossless 
oscillation when the bound state exists; otherwise, the oscillation will decay completely. This result 
is compatible with the coherence-incoherence transition occurring in the small-tunneling limit. Our 
work indicates that the QPT physics in the spin-boson model needs further exploration. 
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I. INTRODUCTION 

The spin-boson model (SBM), describing a two- level 
system (spin) coupled to an infinite collection of har- 
monic oscillators acting as the bosonic reservoir, has 
an extensive relevance to physical systems in con- 
densed matter physics, quantum optics, and quantum 
chemistry. ' An interesting aspect of the SBM is its 
quantum phase transition (QPT) 3 from delocalization 
to localization phases with the increase of the cou- 
pling constant a, as a consequence of the competi- 
tion between tunneling and dissipation induced by the 
reservoir. 4-16 The QPT has been extensively studied 
by many different methods, including a path integral 
method under a non-interacting blip approximation, 4 a 
variational method based on unitary transformation, 5 ' 6 a 
numerical renormalization group, 7 quantum Monte Carlo 
simulation, 8 and direct numerical diagonalization in dif- 
ferent bases, 9,10 A consensus of these methods is that, 
for the ohmic spectrum, the delocalized-localized QPT 
is of Kosterlitz-Thouless type 1,4 and occurs at a = 1 
in the small-tunneling limit. Several experiments have 
been suggested to observe this reservoir-induced QPT; 
these include the use of mesoscopic metal rings, 17 single- 
electron transistors with electromagnetic noise, 18 and 
cold atoms in optical lattices. 19 

While the physics in the localized regime is easy to 
understand, the delocalized regime exhibits many inter- 
esting phenomena. 10-13 For example, it was found that 
in the small-tunneling limit the spin dynamics in the 
delocalized regime shows damped coherent oscillation 
when a < 1/2 and incoherent relaxation for a > 1/2, 
a phenomenon referred to as the coherence-incoherence 
transition. 4 In fact, the understanding of this dynamical 
transition and the physics in the regime of 1/2 < a < 1 
are not clear. Here we focus on the delocalized regime 
and study analytically the ground-state property of the 
SBM by the variational method. Our starting point is to 
examine whether a bound state between the spin and its 
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FIG. 1: (Color online) The phase diagram of the ohmic SBM 
in different parameter regimes, (a) The delocalized-localized 
QPT and (b) the novel QPT found in the delocalized regime, 
which is separated into phases with and without a bound 
state. 



reservoir can form. This idea is inspired by recent exper- 
imental and theoretical developments indicating that the 
formation of a bound state between a two-level system 
and its reservoir can induce decoherence suppression of 
the system in a structured reservoir. 20-23 We find that be- 
sides the well-known delocalized-localized transition, as 
depicted in Fig. 1(a), an alternative QPT is observed in 
the delocalized regime, which accompanies the formation 
of a bound state between the spin and its reservoir, as in- 
dicated in Fig. 1(b). This QPT distinguishes the dynam- 
ics of the system from complete decoherence to decoher- 
ence suppression when the initial state of the reservoir is 
a multimodc coherent state. In the small-tunneling limit 
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this QPT happens at a = 1/2, which is compatible with 
the coherence- incoherence transition point. Figure 1 also 
shows the overall phase diagram for the finite tunneling 
amplitude. Our analytical formulation provides a clear 
physical picture of this QPT and a unified description of 
the QPT in the ohmic SBM. 

Our paper is organized as follows. In Sec. II, we intro- 
duce the SBM and its reduction by the variational ap- 
proach. In Sec. Ill, the novel QPT occurring in the con- 
ventional delocalized-phase regime is revealed. In Sec. IV 
we study the dynamical consequence of this novel QPT 
and reveal an interesting dynamical transition at the crit- 
ical point of the QPT. Finally, discussion and summary 
are given in Sec. V. 



II. THE SBM 



the SBM. 5,6 If the tunneling amplitude is renormalizod 
to zero, the system is in the localized phase and the dy- 
namics is trivial. In contrast, if the rcnomalized tun- 
neling amplitude is nonzero, then the system is in the 
delocalized phase, which displays some interesting dy- 
namical behaviors such as damped coherent oscillation 
and incoherent relaxation. 11 In Fig. 1(a) we reproduce 
the delocalized-localized transition for finite tunneling 
amplitude. Note that Eqs. (2) take into account self- 
consistently the effect of counter-rotating terms in the 
SBM. The unitary transformation separates automati- 
cally the one-excitation and excitation-conservative tran- 
sition H[ from the multiexcitation and nondiagonal tran- 
sition H' 2 . It has been proved that in low-temperature 
and weak-coupling regimes the higher-order perturba- 
tion terms H' 2 can be neglected. 6 Then the transformed 
Hamiltonian has the form H 1 « H + H[ = i? c ff- 



The spin-boson model is defined as 1,4 

H=^ t - + ]T uj k b\b k + £ | (fet + bk)(7z , (1) 
k k 

where e and A are, respectively, the energy difference 
and tunneling amplitude between the two spin states, 
b k (b k ) is the creation (annihilation) operator of the kth 
mode of the reservoir with frequency u>k, and g k is 
the coupling strength between the spin and its reser- 
voir, which is further modeled by the spectral density 
J(uj) = J2 k gl5(ui — Wfe) = 2au\~ s u) s Q{uj c — ui). Here uj c 
is the cutoff frequency, a is the dimcnsionless coupling 
constant, Q(x) is the step function, and s classifies the 
reservoir as subohmic if < s < 1, ohmic if s = 1, and 
superohmic if s > l. 4 In the present work we focus ex- 
plicitly on the ohmic spectrum, which is experimentally 
most relevant, and unbiased case e = 0. 

To study the ground-state property, it is convenient 
to make a unitary transformation H' = UHW where 
U = oq>E fc A fc (&£ ~ b k )* x ]e X p[-ifa y ) and X k = f^. 
Explicitly, H' can be written as H' = H +H[+H 2 , with 
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where C = 4^ i s a constant, a± are the 

raising and lowering operators in the a z basis, x = 
Sfc :2 fr^(^l — bk), v k = ,?A J fc ^ fc acts as a renormalized 
coupling strength, and rj = exp[— t^t] is a renormal- 
ized factor for the tunneling. Following the variational 
method, 5,6,24 the transformation parameters = ^^7-5 
can be determined by minimizing the Bogoliubov-Pcicrls 
free energy. The renormalized factor r\ has been used suc- 
cessfully to characterize the delocalized-localized QPT in 



III. NOVEL QPT IN THE DELOCALIZED 
REGIME 

It is interesting to note that H e s has the same struc- 
ture as the well-known "rotating-wave" Hamiltonian in 
quantum optics, 25 where the total excitation number 
N = J2 k bj.b k + (7 + (7_ is conserved. For this kind of 
Hamiltonian, the Hilbert space is split into the direct sum 
of the subspaces with definite N. In this situation one 
can naively deem that the eigenstate \ipo) = |— , {0^}), a 
tensor product of the respective ground state of the two 
subsystems in zero-excitation subspace with eigenvalue 
Eq = — 5A77, is the ground state of the whole system. 6 
Is this always true? Let us examine the eigen solution of 
H c ff . The eigenstate of i? ff in single-excitation subspace 
can be expanded as \tpi) = c |+,0 fc ) + Y,k c k\-,W, 
where represents the bosonic state with only one 
excitation in the Mh mode. From the eigen equation 
H e s\tpi) = E± one can obtain a transcendental 

equation of the eigenvalue E\ , 
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Ei, (3) 



If Eq. (3) has a real root, one can claim that the system 
exhibits a bound state. From the analytic property of 
y(Ei), one can find that Eq. (3) has one and only one real 
root in the regime (—00, — when y(— ^p) < — is 
fulfilled. For the ohmic spectrum, the condition explicitly 
takes the form 
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where a c is the critical point of forming the bound state. 
It is remarkable that E\ is even smaller than the eigen- 
value Eq. Therefore, the ground-state energy is 



Eq — C, a < a c 
Ei — C, a > a c 



(5) 
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This indicates that, with the formation of bound state, 
the ground state is no longer \(po), but \tpi)- Physically, 
such a sudden change of the ground-state structure man- 
ifests clearly the occurrence of QPT in the system. Wc 
can readily verify that the neglected higher-order pertur- 
bation terms H' 2 give no contribution in the two eigen 
bases, i.e. {<pi\H' 2 \ipj} = — 0,1), which means 
that the neglected term H' 2 has no impact on such level- 
crossing-causcd QPT. This in turn validates our approx- 
imation. We can also verify that the eigenstates of H c ff 
in the subspaces N > 2 actually have larger eigenvalues 
than Ei. This implies that the the higher-boson states 
may not become the ground state. 




FIG. 2: (Color online) Ground-state energy (a) and its first 
derivative (b), both of which are dimensionalized to uj c , as a 
function of coupling constant a. 

As we know, the unitary transformation docs not 
change the eigen spectrum of the system, so the QPT 
occurring in H e g also means a QPT in the SBM. At zero 
temperature, the nonanalyticity of the ground-state en- 
ergy is directly connected to the QPT. The first (or nth) 
order QPT is characterized by the discontinuity in the 
first (or nth) derivative of the ground-state energy To 
further characterize this novel QPT, we plot in Fig. 2 
the ground-state energy and its first derivative as a func- 
tion of a by numerically evaluating Eq. (5). We can see 
that, although E g is continuous, its first derivative suf- 
fers a sudden drop at the critical point a c . This implies 
that the QPT is first-order, which can also be verified 
by the fact that the orthogonality of the ground states 
at two sides of a c will cause the ground-state fidelity to 
suddenly drop to zero. 9 We also plot in Fig. 1(b) the 
overall phase diagram of the SBM in the finite tunneling 
amplitude. With the increase of the bare tunneling am- 
plitude A, the critical point a c increases. In the small-A 
limit, a lower bound a c = | exists. The novel QPT oc- 
curring in the conventional delocalized phase also falls 
in the valid regime of our approximation, i.e. the weak- 
coupling regime. This also validates our approximation. 



IV. DYNAMICAL CONSEQUENCE 

What is the dynamical consequence of this novel QPT? 
Some works have shown that the formation of a bound 
state in a dissipative qubit system under a rotating-wave 
approximation is accompanied by decoherence suppres- 



sion of the qubit. 23,26 This is understandable based on 
the fact that the bound state, as a stationary state of the 
whole system, has a vanishing decay rate and the coher- 
ence contained in it would be preserved during the time 
evolution. In a similar manner we also expect that the 
dynamics in our SBM here is qualitatively changed with 
the occurrence of the novel QPT. In the phase without 
a bound state, the dynamics of the system exhibits com- 
plete decoherence, while in the phase with a bound state 
it exhibits decoherence suppression. 
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FIG. 3: (Color online) Time evolution of Pz(t) under differ- 
ent coupling constants a. The other parameter used here is 
A/u} c = 0.1. The cases of a = 0.05 and 0.25 are in the phase 
without a bound state, and the case of a = 0.55 is in the 
phase with a bound state. 

To confirm our above expectation, we study the 
noncquilibrium dynamics of the system, which is usually 
characterized by the quantity P z (t) = (^(t)|<7 z |^(t)), 
with | >!/(£)) being the time-dependent state of the whole 
system. 1,4 We consider the initial state to be 



i*(o)>=i+>®n 



,-X k (bl-b h ) 



\0k), 



(6) 



which is a product state of the excited state in the a z 
basis of the spin and a multimode coherent state of the 
reservoir. It is noted that such an initial state is dif- 
ferent from the one considered in Refs. 10-13. It is 
straightforward to verify that under the unitary transfor- 
mation U, P z (t) = ($f'(t)\a x \iS>'(t)), where the evolution 
of \i&'(t)) = U\^(t)) is governed by H e s under the ini- 
tial condition |tf'(0)) = \+ x , {0 fc }) with |+ x > = l+> ^~ ) 
being the eigenstate of a x . The evolution of the reduced 
density matrix p' = Tr R (\$>' ){'$>'[) of the spin is governed 
exactly by the following master equation: 25 

P'(t) = - i ^[cT + CT_,p'(t)] + ^[2cT_ P '(t)cT + 



(7) 



where Q(t) = -2Im[f&] + A?/ and j(t) = -2Re[|&]. 
Physically, f2(i) plays the role of time-dependent shifted 
frequency and j(t) that of time-dependent decay rate. 
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The time-dependent factor c(t) satisfies 

c(*) + 4r c (*)+ / c(T)/(*-r)dr = 0, c(0) = l, (8) 
* Jo 

in which the kernel function f(x) = J °° J' {ui)e~ tu ' x duj 
connects to the renormalized spectral density J'(oj) = 
Ylk^k^i 1 ^ ~ By solving Eq. (7) numerically, we 

can get P z {t) = Tr s[p' (t)a x ], as shown unambiguously 
in Fig. 3. We can see that P z (t) decays to zero asymp- 
totically in an oscillatory way when the bound state is 
absent (a = 0.25 or 0.35), while it shows lossless oscil- 
lation (a = 0.55) once the bound state is formed. This 
indicates clearly that the decoherence really can be sup- 
pressed in the phase with the bound state. This dynami- 
cal behavior in turn confirms the existence of this bound- 
state-induced QPT in the SBM. The sudden transition 
from complete decoherence to decoherence suppression 
is a result of QPT occurring in the model itself. It is 
the abrupt change of the ground-state structure near the 
critical point of the QPT that induces the qualitative 
difference in the dynamics of the SBM. 

On the other hand, a widely studied case is the 
noncquilibrium dynamics of the spin system when the 
reservoir is initially in a vacuum state. 10-13 It was shown 
that the dynamics changes from damped coherent os- 
cillation to incoherent relaxation, i.e. the so-called 
coherence-incoherence transition, at the point a = i in 
the small-A limit. 4 As a dynamical behavior, this is of 
course not a generic phenomenon for any initial reser- 
voir state. Therefore, the coherence-incoherence transi- 
tion itself does not suffice to deduce the occurrence of the 
intrinsic QPT of the model. Then what is the physical 
reason for this dynamical transition? We notice that the 
transition point matches perfectly well with the critical 
point ct c in the QPT in the small-A limit. Thus, it is 
reasonable to conjecture that the coherence-incoherence 
transition is also a dynamical consequence of this novel 
QPT at the critical point a c . 

The different dynamical consequences are actually due 
to the different reservoir properties. To the initial 
vacuum state of the reservoir, the QPT induces the 
coherence-incoherence dynamical transition, while to the 
multimodc coherent state, it induces the dynamical tran- 
sition from oscillatory damping to lossless oscillation. 
The distinctive character of the coherent state Eq. (6) 
resides in that it is a vacuum state for the reservoir part 
in the H c h representation so that only the two lower eigen 
subspaccs with N = and 1 of the whole system are in- 
volved in the dynamics. This means that only one of an 
arbitrary infinite number of modes of the reservoir with 
one boson affects the interaction with the spin. As a 
result the reservoir shows more coherent interplay with 
the spin and efficiently prolongs the coherence time of 
the spin, especially when the bound state is formed. On 
the other hand, if the reservoir is initially in the vacuum 



state, i.e., |W(0)) = |+i{0fe}), the dynamical behavior of 
lossless oscillation can no longer appear. This is because, 
in this case, the reservoir in H e Q representation corre- 
sponds to Yik l-^fc) with the total boson number being 



E 

k 



J(w) 



4(cj + Ar])' 



-duj, 



which implies that a process of a large number of bosons 
is involved in the dynamics. Such a large number of 
bosons distributed among the infinite reservoir modes 
inevitably induces a completely out-of-phase interaction 
with the spin. Consequently, the lossless oscillation can- 
not occur for the initial vacuum state even in the phase 
with the bound state. In the phase without the bound 
state, Afe and n are not so large because of the small a in 
this phase. Then our results coincide qualitatively with 
the decaying oscillation in Ref. 10. 



V. DISCUSSION AND SUMMARY 

In conclusion, using the variational method, we have 
found a novel QPT in the delocalizcd regime of the SBM. 
This QPT separates the delocalized regime into phases 
with and without a bound state. In these two phases, 
the decoherence dynamics of the spin shows qualitative 
differences. In particular, if the initial reservoir is in the 
multimode coherent state, the spin dynamics in the phase 
with a bound state exhibits a lossless oscillation. In con- 
trast, in the phase without a bound state, the oscilla- 
tion will decay completely to zero. Our purely analyti- 
cal treatment provides a unified microscopic description 
of the QPT and renders a helpful understanding of the 
rich physics in the SBM. Further extension of the present 
work to other spectra is straightforward. As a final re- 
mark, the dynamics discussed in the present work only in- 
volves the zero-temperature case. It is also interesting to 
discuss the dynamical consequence caused by this QPT 
when the reservoir is initially in thermal equilibrium. 27-29 
However, in this case it is very difficult to obtain the ex- 
act dynamics because the initial state involves the eigen 
bases of all the subspaces with definite excitation number 
from zero to infinity. Therefore, some approximations are 
necessary. This problem deserves further exploration in 
the future. 
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